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ABSTRACT 



Context. Paczynski realized that a properly chosen gravitational potential may accurately model (in a "pseudo Newtonian" theory) 
general relativistic effects that determine motion of matter near a non-rotating black hole. Paczyhski's choice, known today as the 
"Paczyhski-Wiita potential", proved to be very practical. It was used by numerous researchers in the black hole accretion theory, and 
became a standard tool in relativistic astrophysics. 

Aims. The model is an example of Paczyhski's admired ability to invent "out of nowhere" simple ideas that were brilliant, deep 
and useful. Paczynski has guessed intuitively the form of the potential. However, it could be also derived by a a step-by-step formal 
procedure. I show the derivation here. 

Methods. My derivation is based on a standard definition of the relativistic "effective potential" in the Schwarzschild spacetime. 
Results. The relativistic effective potential may be uniquely divided into its "gravitational" and "centrifugal" part. The gravitational 
part differs from the Paczyhski-Wiita potential only by a constant. 
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1 . A brief history of the potential 

|Abramowicz, Jaroszyriski & Sikora 1978| working in Paczyii- 
ski's research group in Warsaw, found a practical mathemati- 
cal scheme to construct fully relativistic models of thick accre- 
tion disks, known today as "Polish doughnuts". The models dis- 
played several astrophysically interesting features, among them 
seriously super-Eddington luminosities, long and narrow fun- 
nels that could collimate jets, and a self-crossing "Roche lobe" 
equipotential. The scheme developed in Warsaw was suitable for 
both analytic and numerical calculations. 

Some properties of the thick disks were obviously connected 
to strong-field effects of general relativity. Our leader Bohdan 
Paczynski, who was not familiar with the technicalities of gen- 
eral relativity, asked me to find a Newtonian way to describe 
these effects. I was rather unhappy about Paczyiiski's request, 
because initially I wrongly imagined that the only method ade- 
quate for the task should be the post-Newtonian scheme. It uses 
tedious, long and boring expansions. I was desperately working, 
producing longer and longer formulae, when one day Paczynski 
came to my office, and said "Stop working on that. I found the 
solution." And he showed me his solution — a NewtoniarQ po- 
tential. 



the two most important radii that characterize circular Keplerian 
orbits, the radius of marginally stable orbit r,„s (i.e. ISCO) and 
the radius of marginally bound orbit r„,b have the same values in 
Newton's gravity with his potential ^ as in Einstein's gravity in 
the Schwarzschild metric. 
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where r is the spherical radius, M is the mass of the black hole, 
and rc is its gravitational radius. Paczynski had checked that 



' It is often called "pseudo Newtonian" to stress that it does not obey 
the Poisson equation. However, when the external gravity is fixed, (as it 
is in the black hole accretion theory) "pseudo Newtonian" is practically 
equivalent to "Newtonian". This is why I am using both terms here. 



This was a brilliant display of the qualities of Paczyhski's mind: 
he just guessed the right, simple and powerful solution to the 
problem. His solution immediately proved to be very practical. 
Shortly afterwards, )Paczyriski & Wiita 1980 used ([TJ to numeri- 
cally calculate models of thick disks. The models differed from 
these calculated with the full strength of general relativity by 
only a few percent. Later, this opened a flood gate: numerous 
authors used the Paczyhski-Wiita potential in their calculations 
of black hole accretion flows. The potential is so remarkably suc- 
cessful, that some researchers use it even outside its obvious lim- 
its of applicability: (a) for rotating black holes; and this is wrong 
because ([TJ does not include the Lense-ThiiTing effect, (b) for 
self-gravitating fluids; and this is wrong because V^O it 0. 



2. A "derivation" of the Paczyhski-Wiita potential 

Why is the Newtonian Paczyhski-Wiita potential ([T) such an ac- 
curate model of the strong relativistic effects? Should this be 
considered a fortunate, unexpected coincidence, or could one 
"derive" the potential from the first principles of Einstein's gen- 
eral relativity? I remember discussing this question shortly with 
Thibault Damour in the late 1970. Although we were convinced 
that the "effective potential approach" should provide such a 
derivation, we have not completed relevant calculations. I sum- 
marize them here. 
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In Newtonian theory, let E be energy, L angular momen- 
tum, 0(r) gravitational potential, and V radial velocity. The or- 
bital motion is often described in terms of the effective potential 
U{r,L)^<b(r) + L^I2r^, 

1 



3. Other potentials, special relativity 



-V' ^E-U(r,L), 
Circular orbits locate at the effective potential extrema, 

or in terms of the gravitational potential O, 
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Let us consider almost circular motion of particles on the 6 = 
Jill "equatorial" plane in the Schwarzschild spacetime. From 
the radial component of the four-velocity u'' let us construct a 
positive small quantity = (Ur)^ g'''' «: 1. It is known that 
Mf and M0 are constants of motion, therefore L = -u^/u, is 
also a constant of motion (the specific angular momentum). The 
1 = M/ Uk g''' - (ut)^g" + {Utpfg'^'^ + {urfg'''' condition may be 
written in the form. 



^ ln( 1 + y2) = ln(M,) + ^ In [g" + g^^\ (6) 

Expansion of the left-hand side yields 12. One also defines 
E = ln(M,) and, 

Uir,L)^-]^\n[g" +L^g^^\. (7) 

This brings equation (|6]l into a form identical with the Newtonian 
formula (O. Thus, the Newtonian condition (HI for the vanish- 
ing derivative of the effective potential may be applied to the 
relativistic effective potential which gives 
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Because at the equatorial plane g'^ 
this may be written in the form, 

dl GM\ 
dr\ r - rc] 

Comparing Newton's equation (|5]l with Einstein's equation (|9]l, 
we see that the gravitational potential in both equations has to 
have the same Paczynski-Wiita form (HJ. In deriving equation 
(|9]l we used, 

= — ^ = -H 1 =2(D(r)H- 1. (10) 
r- Kg r-rc 

Thus, the Keplerian angular momentum derived (in the 
Schwarzschild spacetime) according to Einstein's theory, and 
derived with the Paczynski-Wiita potential are both given by the 
same formula, 

^2 GMr^ ni) 

In Newton's theory the angular momentum L and angular veloc- 
ity Q are connected by L = Qr^, but in Schwarzschild geometry 
by, L - Qr^/(1 - rdr). Therefore, the Keplerian angular ve- 
locity calculated in Schwarzschild geometry and in Paczynski- 
Wiita potential are not the same, 
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Nowak & Wagoner 1991 [found that the potential given by a fit- 



ting formula 
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reproduces the angular velocity Q(r) and the radial epicyclic fre- 
quency K(r) better than the Paczyfiski-Wiita potential. The fit- 
ting formula used bv lKluzniak & Lee 2002l reproduces the ratio 
K{r)IQ.{r) of these frequencies exactly: 

ISemerak & Karas 1 999 discussed the Newtonian potential suit- 
able for modeling the gravity of the Kerr black hole, including 
the Lense-Thirring effect. 

Neither any of these three potentials, nor a few other po- 
tential introduced by some other authors, became popular 
Nowadays much more astrophysicists know Einstein's general 
relativity than in the late 1970s, but quotations to Paczynski- 
Wiita potential show no sign of declinfl 

Velocities of matter calculated with the Paczyriski-Wiita po- 
tential could exceed the light speed. This creates a serious prob- 
lem when one calculates the observed appearance of matter (e.g. 
spectra) by the method of ray tracing. I Abramowicz & al. 19961 
found a fix for this problem by showing how to incorporate the 
effects of special relativity into the Paczynski-Wiita scheme: one 
should interpret the "true" physical velocities in terms of the cal- 
culated ones by Vcai - V,„,/(l - V,^„/c^). Here V(,„) denotes each 
of the three components of the velocity, i.e. V,-, Vg, V^. 

4. Conclusions 

The Paczyriski-Wiita potential ([T]i accurately models in the 
Newtonian theory general relativistic effects that determine mo- 
tion of matter near a non-rotating black hole. 

The Paczyriski-Wiita potential is neither an approximation of 
relativistic gravity, nor a fitting formula. Instead, it is a unique 
("pseudo") Newtonian model of the gravity of a non-rotating 
black hole. It reproduces exactly: 

- The location of the marginally stable orbit r„,j = ISCO. 

- The location of the marginally bound orbit rmt,- 

- The form of the Keplerian angular momentum L{r). 

It also reproduces accurately, but not exactly, the form of the 
Keplerian angular velocity Q.(r), and of the radial epicyclic fre- 
quency K{r). 
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